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Abstract 

I , We study rotating boson stars in five spacetime dimensions. The boson fields consist of a 

^Jy complex doublet scalar field. Considering boson stars rotating in two orthogonal planes with 

both angular momenta of equal magnitude, a special ansatz for the boson field and the metric 
allows for solutions with nontrivial dependence on the radial coordinate only. The charge of 
the scalar field equals the sum of the angular momenta. The rotating boson stars are globally 
regular and asymptotically flat. For our choice of a sixtic potential the rotating boson star 
solutions possess a flat spacetime limit. We study the solutions in flat and curved spacetime. 
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Boson stars are globally regular configurations of self-gravitating boson fields 1, 2, 3,4.. When the 
complex scalar field has a suitable self- interaction, the boson stars possess a flat spacetime limit, 
representing non-topological solitons or Q-balls, OE]. These are stationary localized solutions with 
a finite mass, whose charge Q corresponds to their particle number. Supersymmetric extensions 
of the Standard Model also possess Q-ball solutions [7|. In this case several scalar fields interact 
via complicated potentials. It was shown that cubic interaction terms that result from Yukawa 
couplings in the superpotential and supersymmetry breaking terms lead to the existence of Q-balls 
with non- vanishing baryon or lepton number or electric charge. These supersymmetric Q-balls have 
been considered recently as possible candidates for baryonic dark matter [5] and their astrophysical 
implications have been discussed [9]. Supersymmetric Q-balls have been constructed numerically 
in [10]. 

The complex scalar field of boson stars and Q-balls has a harmonic time-dependence with 
frequency u) s . In 4-dimensional flat spacetime, the allowed frequency range for Q-balls w m i n < 
lo s < w m ax is solely determined by the properties of the potential describing the self-interaction 
[6] . This frequency range holds for both non-rotating and rotating Q-balls [11] [TSl [13] . For 
4-dimensional boson stars the maximal frequency w max remains the same as in flat spacetime, 
whereas the minimal frequency depends on the strength of the gravitational coupling constant 
and increases with increasing gravitational coupling. 
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For Q-balls in 4 dimensions, the mass and charge diverge at both ends of the frequency interval 
and assume their minimal value at a critical frequency cj cr , where the classical stability of the 
solutions changes [SJ|B]. For uj s < uj cr , i.e., below the critical value, the solutions are classically 
stable. As long as their mass is smaller than the mass of Q free bosons they are also quantum 
mechanically stable 5,- In constrast, above w cr the solutions are unstable. 

The coupling to gravity modifies and enriches this simple pattern. Thus when w s — ¥ w max , the 
mass and charge of boson stars tend to zero, when the spacetime is 4-dimensional. Moreover, in 
the lower frequency range, the boson star solutions are not uniquely determined by their frequency. 
Instead a spiral-like frequency dependence of the charge and the mass is observed, where the charge 
and mass approach finite limiting values at the center of the spiral. Thus a plethora of boson star 
branches appears, typically featuring a single classically stable branch |14) . 

Remarkably, for charged boson fields with a F-shaped scalar potential, also boson shells can 
arise [HI EI]. These may harbour black holes, violating black hole uniqueness [17] . 

The presence of rotation leads to an interesting feature of these solutions, namely their angular 
momentum is quantized and corresponds to an integer N times the charge, J = NQ 118] . 
Rotating Q-balls and boson stars with N = 1 thus have the smallest angular momentum J for a 
given charge Q. While non- rotating Q-balls and boson stars are cohomogeneity-1, the presence 
of rotation includes an angular dependence and thus leads to a cohomogeneity-2 problem in 4 
dimensions. 

With the advent of string theory the exploration of higher dimensions has yielded many new 
insights for localized objects and their properties. For black holes, in particular, it became clear, 
that many features of black holes in 4 dimensions do not extend to higher dimensions [TJj]. It 
therefore appears eligible, to consider also other localized objects in higher dimensions, such as 
boson stars, which still represent alternatives to black holes [3D]. A first effort in this direction was 
undertaken by Astefanesei and Radu |21] , who addressed non- rotating boson stars with a negative 
cosmological constant in 5 dimensions, and by Prikas |22] who studied electrically charged non- 
rotating boson stars in various dimensions in asymptotically flat as well as Anti de-Sitter spacetime. 
Here we consider rotating asymptotically flat boson stars and their flat spacetime counterparts. 

In general, rotating objects in D dimensions possess [(D— l)/2] independent angular momenta 
[2 31 . Thus in 5 dimensions, there are two independent angular momenta associated with the two 
orthogonal planes of rotation. Generic rotating solutions in 5 dimensions then depend on two 
coordinates, just as in 4 dimensions. However, when both angular momenta are chosen equal, 
J\ = J2 = J, the symmetry of the solutions may be enhanced. For black holes one then obtains a 
cohomogcneity-1 problem, with dependence only on the radial coordinate, also in the presence of 
Abelian gauge fields and dilaton fields |24) . 

Here we show, that with an appropriate Ansatz for the scalar fields, the problem can be reduced 
to cohomogeneity-1 also for rotating Q-balls and boson stars in 5 dimensions, when both angular 
momenta are chosen equal. For that purpose we introduce a complex doublet of scalar fields as 
present, for instance, in the electroweak sector of the standard model. Moreover, we show that 
the quantization of the angular momentum also carries over, yielding for the lowest such rotating 
solutions the relation 2\J\ = Q. 

The paper is organized as follows. In section 2 we present the action, the ansatz and the 
expressions for the physical properties of the rotating Q-balls and boson stars in 5 dimensions. 
In section 3 we exhibit the numerical results for the solutions in flat spacetime, comparing non- 
rotating and rotating Q-ball solutions, and discussing their stability. The effects of the coupling 
to gravity is the topic of section 4, where we present the non-rotating and rotating boson star 
solutions. We give our conclusions in section 5. 
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2 Action, Ansatz and Charges 
2.1 Action 

In 4 dimensions, Q-balls and boson stars with 2 complex scalar fields have been considered before, 
leading to interesting phenomena due to their interaction |25j . Here we introduce a complex 
doublet scalar held $, in order to obtain rotating cohomogeneity-1 Q-balls and boson stars in 5 
dimensions. 

Thus we consider the following action of the self-interacting complex doublet scalar held $ 
coupled minimally to Einstein gravity in 5 spacetime dimensions 



S = 



16ttG 



^gd 5 x, (i) 



with curvature scalar R and Newton's constant (in 5 dimensions) G and self-interaction potential 
U . Here t denotes the complex transpose, and |$| 2 = 
The scalar potential U is chosen as 

{/(|$|)=A(|$| 6 -a|<I>| 4 + &|<i>| 2 ), (2) 

where A, a, b are constants. The choice of potential allows for the existence of nontopological 
soliton solutions, i.e., Q-balls, in the absence of gravity. This self-interaction of the scalar field 
has an attractive component, and the potential has a minimum, U(0) — 0, at $ = and a second 
minimum at some finite value of |$|. The boson mass is determined by the quadratic term and 
given by tob = n/A6. 

Variation of the action with respect to the metric and the scalar fields leads, respectively, to 
the Einstein equations 

G M „ = Rfj, u — -g^R — 87rGT M „ (3) 

with stress-energy tensor 

. dL M 



dg» v 



= -\g„ v ((8 a $)*(dp$) + (dp^(d a $)) g af} + (0 p *)t(a,*) + (d^)Hd^) 

-\g^U(m) 7 (4) 
and the matter field equation, 

-±=B, (V- 9 0^) = ^L*. (5) 

2.2 Ansatz 

To construct stationary rotating solutions with two equal angular momenta we employ bi-azimuthal 
isotropic coordinates [2SI HI] 

ds 2 = g^dx^dx" = -fdt 2 + — (dr 2 + r 2 d9 2 ) + —j^r 2 sin2 6 cos2 9 ( d( Pi ~ d W 

+ -r 2 ( sin 2 etdcf! - -dt) 2 + cos 2 6(d(p 2 - -dt) 2 ) . (6) 
j \ r r ) 

This metric has 3 commuting Killing vector fields 

i = d t , r]i=d tpi , f] 2 = <9 V2 , (7) 

where £ denotes an asymptotically timelike Killing vector field and r\i denote asymptotically space- 
like Killing vector fields. The 4 metric functions /, m, n and uj are functions of the radial variable r 
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only. In the non-rotating case, the metric simplifies, since m = n and u = 0. In the flat spacctimc 
limit, all metric functions become trivial, f = m = n= l,u) = 0. 

For the complex doublet boson field we assume a harmonic time-dependence with frequency 
lo s . The Ansatz then has the form 

$ = cj>(r)e iUet ^, (8) 

where $ denotes a complex unit two vector. To specify this unit two vector <l, we distinguish 
between non-rotating and rotating solutions. For non-rotating solutions we choose 



<!>=( J ] (!)) 



whereas for rotating solutions we choose 

<i>=f sin ^M (10) 

This choice ensures single- valuedness of the scalar field, since 

$(v?i,¥>2) =$(2?r + yi,27T + <p 2 ) • (11) 

The specific angular dependence yields a cohomogeneity-1 system. 

The choice of appropriate boundary conditions must guarantee that the solutions are globally 
regular and asymptotically flat, and that they possess a finite energy and finite energy density. 
We therefore impose at infinity the boundary conditions 

f\r— oo — ^|r— oo — ^|r— oo 1 : ^|r— 00 0|r— 00 — > (1^) 

thus the scalar field assumes its vacuum value $ = 0. At the origin, on the other hand, regularity 
requires the boundary conditions 



together with 



<9 r /| r=0 = d r m\ r=0 = (13) 
d r 0\ r =o = (14) 



for the non-rotating solutions, and 

d r n\ r=0 = , oj\ r =o = 4>\r=o = (15) 

for the rotating solutions. 
2.3 Charges 

The mass M and the angular momenta Ji, i — 1,2, can be obtained from their respective Komar 
expressions 



and 



\ I R^nOfdV , (16) 



E 



Here £ denotes an asymptotically flat spacelike hypersurface, is normal to £ with n^n^ = — 1, 
dV is the natural volume element on E, and £ and rji denote the Killing vector fields. Replacing 
the Ricci tensor by the stress-energy tensor yields 

M = \ jT (t„„ - \g^T^j n^CdV , (18) 
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and 

Ji = ~ (t»u - \g^T^j n»rffdV . (19) 

The mass M and the equal angular momenta J = J\ = J2 can also be obtained directly from 
the asymptotic expansion for the metric 



8 GM 4GJ „{1 



/ = 1 -CT + 0(- I ), (20) 



The conserved charge Q is associated with the scalar field $, since the Lagrangian density is 
invariant under the global phase transformation $ — > $e la . This leads to the conserved current 

f = -i - 0"$t$) , f.^ = , (21) 

with the conserved charge Q 

Q = - J .fV=gd 4 x. (22) 
Inserting the Ansatz into the expressions for the angular momenta Ji 

Ji = -2tt 2 J v/W"^ (w. + ^) 4> 2 r Z dr (23) 

and the conserved charge Q 

Q = 4ir 2 J v/W^ (u. + 5?) 4> 2 r 3 dr (24) 
we obtain the quantization condition 

J2\Ji\=2\J\=Q . (25) 

i 

3 Solutions in flat spacetime: Q-balls 

In flat spacetime, i.e. for G = the model possesses Q-ball solutions. In this limit the metric 
functions are trivial f = m = n = 1, ui = 0. 

3.1 Q-ball properties 

The mass M of non-rotating Q-balls is given by 



M = 2vr 2 / T tt r 3 dr = 2tt 2 / [co 2 4> 2 + 4>' 2 + U(<f>)] r 3 dr , (26) 
Jo Jo 

where the prime denotes differentiation with respect to r, and their charge Q is given by 

Q(uj s ) = 4ir 2 u s <p 2 r 3 dr. (27) 
Jo 

As in 4 dimensions [11J . Q-ball solutions exist in the frequency range 

w min <uj2 s < W^ ax , (28) 

where 

«*L* = ^"(0) = Xb = m%, (29) 

and 

< in = mm [U(t)/0 2 ] =\(b-?p) . (30) 
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Thus the maximal value is determined only by the boson mass tub- 

For rotating Q-balls the mass receives an additional contribution and becomes 



M = 2tt 2 / u 2 cj> 2 + 3-V + (f>' 2 + U(6) 
Jo L r 

The charge and the frequency range remain unchanged. 



r 3 dr , (31) 



3.2 Numerical results 

For the numerical calculations we introduce the compactified radial coordinate r = r/(l + r) [27] . 
We employ a collocation method for boundary-value ordinary differential equations, equipped 
with an adaptive mesh selection procedure [28]. Typical mesh sizes include 10 2 — 10 3 points. The 
solutions have a relative accuracy of 10 -6 . 

In the calculations of the 5-dimensional solutions we choose the same set of potential parameters 
as employed previously in 4 dimensions [TTJ [T2J [T3] 

A=l, a = 2, 6 = 1.1. (32) 

This yields for the boson mass me = s/hl. We denote the gravitational coupling constant by 
k = 8nG. 

In 4 dimensions, when fixing the value of lo s in the allowed range, one obtains a sequence of 
Q-ball solutions, consisting of the fundamental Q-ball and its radial excitations [TT]. The boson 
function of the fundamental Q-ball has no nodes, while it has k nodes for the fc-th radial excitation. 
Here we focus only on the fundamental solutions for Q-balls in 5 dimensions. 

We exhibit the mass M and charge Q of non-rotating and rotating Q-balls in Fig. 1(a) versus 
the frequency u s . At the limits of the interval, i.e., for uj s — > uj m - m and uj s — > w max , both mass 
and charge diverge. Inbetween, mass and charge assume their minimal values, from where they 
rise monotonically to both sides. For a given frequency u„ the rotating configurations possess 
higher mass and charge. The dependence of the global charges on the frequency thus is completely 
analogous for these 5-dimensional Q-balls as for their 4-dimensional counterparts [12j [13] . 



In Fig. 1(b) properties of the scalar field are exhibited. For non-rotating Q-balls the value 
of the scalar field function at the origin 0(0) is shown versus the frequency, whereas for rotating 
Q-balls the value of the derivative of the scalar field function at the origin d r 4>\ r =o is shown. 
Both tend to zero at the upper frequency limit, while they assume finite values at the lower limit. 
Inbetween they assume a maximal value. 

3.3 Stability 

To address the stability of these solutions, we follow the arguments presented in ref. |TJ. The 



stability of the non-rotating solutions can be read off from Fig. 1(c) where the mass M of the 
Q-balls is exhibited versus their charge Q, together with the mass Mf = tiibQ of Q free bosons. 

We thus observe two branches of Q-balls solutions, a lower branch and an upper branch. On 
the lower branch the solutions are classically stable. As long as M < Mf, i.e., for most of this 
lower branch, they are also quantum mechanically stable. Shortly after crossing the free mass 
Mf, however, the lower branch ends at the critical point, where mass and charge assume their 
minimal values. The upper branch starts from this critical point, and approaches the free mass 
Mf asymptotically. At the critical point the classical stability changes, and the upper branch has 
an unstable mode [T] [Tl] . 



The situation is analogous for rotating Q-balls, as seen in Fig. 1(d) Here the critical point, 
where the classical stability changes, is encountered for larger values of the mass and charge. 
Otherwise, the presence of rotation does not seem to be relevant for the issue of stability. 
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Figure 1: Properties of Q-ball solutions are shown: (a) Mass M and charge Q versus the frequency 
cj s ; the vertical lines correspond to the minimal and maximal values of uj s . (b) Value of the scalar 
field function at the origin <fio = <fi(0) for non-rotating Q-balls, and the value of the derivative of 
the scalar field function at the origin <fii = d r 4>\ r= o for rotating Q-balls. The asterisks indicate the 
extrapolated limits for u> s — > uj m j n . (c) Mass M versus the charge Q for non-rotating Q-balls; also 
shown is the mass of Q free bosons, (d) Mass M versus the charge Q for rotating Q-balls; also 
shown is the mass of Q free bosons. 
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4 Solutions in curved spacetime: boson stars 

We now consider the effect of gravity on these solutions. The resulting boson stars have many 
features in common with their 4-dimensional counterparts. But in 5 dimensions also new features 
arise. In the following we first discuss the properties of non-rotating boson stars and then turn to 
rotating boson stars. 

4.1 Non-rotating boson stars 
4.1.1 Properties 

For the boson star solutions in 5 dimensions the frequency dependence is in many respects analo- 
gous to the 4-dimensional case. The frequency u> s is also bounded from above by w max , Eq. (|29|) . 
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Figure 2: Properties of the non-rotating boson stars: (a) Mass M versus rotation frequency w s 
for several values of the gravitational coupling constant k. The vertical lines correspond to the 
flat spacetime minimal and universal maximal values of lu s . (b) Same as (a) for the charge Q. 
(c) Same as (a) for the scalar field function <fio = (f)(0) at the origin. The asterisk indicates the 
extrapolated limit uj s — > w m j n for k = 0. (d) Same as (a) for the scaled mass kM. 



because the scalar field exhibits an exponential fall-off asymptotically only for lu s < U) max . More- 
over, the frequency is also no longer bounded by the flat spacetime minimal value ui m i n , Eq. (l30|) . 
because the presence of gravity fundamentally changes the dependence of the mass and charge on 
the frequency, and a spiralling pattern arises at the lower frequency end. 

The frequency dependence of the mass M and charge Q is exhibited in Figs 



2(a) and 2(b) 



for non-rotating boson stars and a sequence of values of the gravitational coupling k. While the 
presence of a spiral is a genuine property of boson stars, the location and the size of the spiral 
depend on the gravitational coupling strength. For very small k the spirals have a small extent in 
Ld s and are located beyond w m j n , thus at smaller values of the frequency than available in the flat 
spacetime limit. With increasing k the spirals then move to larger frequency values and at the 
same time cover a larger range of frequencies. 

We exhibit the frequency dependence of the value of the scalar field at the origin 0(0) of these 
boson star solutions in Fig 2(c) The upper endpoints of (f)(0) in the figure correspond to values 
close to the centers of the spirals, where the numerical procedure was stopped. Thus these upper 
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endpoints do not signify any physical relevance. This is in contrast to the universal lower endpoint 
of 0(0), 0(0) -> for oj s -> w max . 

Let us then inspect the behaviour of the boson star solutions in the limit lo s — > oj max more 
closely. Interestingly, in this limit the mass and the charge of the 5-dimensional boson stars remain 
finite. In 4 dimensions this is not the case, but the mass and charge tend to zero in this limit. 
To obtain further insight into this at first unexpected behaviour, we consider the scaled mass kM 
and scaled charge kQ. The scaled mass is exhibited in Fig. |2(d)| and we note that it assumes a 
universal value for all boson star solutions in the limit lu s —> w max , and so does the scaled charge. 

We can understand this behaviour, when we realize that the metric becomes Minkowskian in 
this limit, thus / = m = 1, while the boson field spreads and tends to zero. At the same time 
the scaled boson field function 

0(f) = 0(f )/0o , with r=^ Q ^y r ,, (33) 

tends in this limit uj s — > w max to a universal function for all k. 

We now express the charge in terms of the scaled boson field function and the scaled coor- 
dinate f 

1 



3C 



Q = 4tt 2 w s - / 2 f 3 df , (34) 

and likewise the mass 



,, 2ir 2 r°° 
M = — 



f 3 df . (35) 



In the limit uj s — > cJ max the value of the scalar field at the origin tends to zero, 0o — > 0, and thus 
the mass tends to 

A 2 poo 

M = — < ax / 2 f 3 df . (36) 

K 

Since the charge tends to 



Q = 4tt 2 uw~ / 2 f 3 df , (37) 
K Jo 

we obtain the limiting relation between mass and charge 

M = u max Q = m B Q ■ (38) 

The explicit limiting values of kM and nQ are determined by the limiting universal function 0. 

To determine 0, we first observe that for all k the ratio (w 2 lax — oj 2 )/0ok5 assumes the same 
value w 2 , say. We then make a perturbative expansion for the metric functions 

/ = e 5v ^ , m = e^ (f) , (39) 

where S is small. From the numerical solutions we find S = 4>qk^ . Substitution in the boson field 
equation and the Einstein equations yields to lowest order in 0o 

■ + »L) . ^ (f 3 ^) = !^ 3 <^iax , A* = v /i ■ (40) 




The solution to this set of equations is exhibited in Fig. 3(a) 

Let us now compare the limit u> s — > a; max for the boson star solutions in 4 and 5 spacetime 
dimensions. In constrast to 5 dimensions, where charge and mass assume finite values as lo s tends 
to its maximal value, they vanish in this limit in 4 dimensions. Analysing the 4-dimensional 
solutions we observe that the scaling behaviour of the boson field function is the same as in 5 
dimensions, see Eq. (|33D . However, the volume integral of the mass and charge in 4 dimensions 



9 



10 



10 



(a) 



(b) 



Figure 3: (a) The function <j>(f) and v(f) are shown for the non-rotating boson stars, (b) The 
same for the rotating boson stars. Also shown is the function w(r). 



contains a power of r less than the volume integral in 5 dimensions. Consequently the powers of 
(j)o do not cancel, 



) (d=4) flm , ^0 

> y ' = »7rw max — 



4> 2 f 2 df 

and the charge vanishes as 4>o tends to zero, i.e., in the limit u> s ■ 



(41) 



4.1.2 Stability 

Since the mass M and the charge Q have an analogous frequency dependence, we obtain an 
intricate cusp structure when we consider the mass as a function of the charge as illustrated in 
Fig. |U As long as gravity is weakly coupled, e.g., k < 0.2, the single cusp present in fiat spacetime 
(associated with the minimal value of the mass and charge) is retained. Here the classically stable 
branch of boson star solutions ends, and a classically unstable branch arises. 

The boson star solutions on this unstable branch have a mass that is larger than the mass of 
Q free bosons. The branch has a finite extent and ends when w s — > w max , where M = Mf = TObQ- 
This behaviour is different from 4 dimensions, where an additional cusp occurs close to u; max , 
because mass and charge tend to zero in the limit. As seen in Fig. 2J the relative extent of this 
classically unstable branch as compared to the classically stable branch increases as the coupling 
k increases. Being barely noticable for k = 0.01, Fig. |4(a) the classically unstable branch is the 
dominant branch (qua extent), when k = 0.2, Fig. 4(d) 



The inlets in Fig. [4] demonstrate the cusp structure due to the presence of the spiral. Whenever 
mass and charge assume extremal values, another cusp arises. According to catastrophe theory, 
however, at each new cusp another unstable mode arises [14] . Thus the boson stars become 
increasingly unstable as the central point of the spiral is approached. 

Turning finally to large values of the gravitational coupling, we observe, that the classically 
stable branch has disappeared completely, since mass and charge no longer exhibit a minimum 
outside the spiral. Moreover, all boson star solutions now have masses above the mass of Q free 
bosons. Thus in this range of k only unstable solutions remain. 

Since for large k gravity dominates, the boson self-interaction terms become negligible, as 
demonstrated explicitly in 4 dimensions [T2]. These solutions then correspond to boson stars 
with only a mass term. In 5 dimensions such solutions were investigated in the presence of a 
negative cosmological constant [21J . Interestingly, the cosmological constant again allows for a 
stable branch, as explicitly demonstrated by Astefanesei and Radu [21] . 
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Figure 4: Mass M versus charge Q for non-rotating boson stars for several values of the gravita- 
tional coupling constant: (a) k — 0.01, (b) k — 0.05, (c) K = 0.1, (d) n — 0.2, (e) k = 0.5, (f) 
k = 1.0; also shown is the mass of Q free bosons. 
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4.2 Rotating boson stars 

Let us now turn to rotating boson stars in 5 dimensions with equal angular momenta. Interestingly, 
their properties are very similar to those of the non-rotating boson stars. 

4.2.1 Properties 
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Figure 5: Properties of rotating boson stars: (a) Mass M versus rotation frequency uj s for several 
values of the gravitational coupling constant k. The vertical lines correspond to the flat spacetime 
minimal and universal maximal values of tu s . (b) Same as (a) for the charge Q. (c) Same as (a) 
for the derivative of the scalar field function <f>\ = d r <j)\ r= o at the origin. The asterisk indicates 
the extrapolated limit u s — > Li m i n for k — 0. (d) Same as (a) for the scaled mass kM. 



The frequency dependence of the mass M and charge Q of rotating boson stars is in general 



rather similar to the case of non-rotating boson stars, as seen in Figs. 5(a) and 5(b) for an analogous 
sequence of values of the gravitational coupling k. The only striking effect of the rotation is that the 
spirals elongate substantially for the somewhat larger values of the coupling strength. However, 
this effect is not unexpected, since we observed it previously also for rotating boson stars in 4 
dimensions [TJl H3] (representing cohomogeneity-2 solutions). 

The frequency dependence of the derivative of the scalar field function at the origin <f>x = 
d r cf>\ r= o, exhibited in Fig. 5(c) nicely illustrates these large spirals by the large oscillations of uj s 
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with increasing 4>i ■ 

Let us now again inspect the behaviour of the boson star solutions in the limit cj s — > w ma x- 
As in the non-rotating case, the mass and the charge of the 5-dimensional rotating boson stars 
remain finite in this limit. The scaled mass kM is exhibited in Fig. 5(d)| and we observe again 



a universal value for the scaled mass and also for the scaled charge for all boson star solutions in 
the limit uj s — > u; max , with the same relation M — m^Q, Eq. (|38p . However, the explicit values of 
kM and kQ differ from the non-rotating case. 

To see where this difference comes from, we follow the same set of arguments as in the non- 
rotating case and introduce a scaled boson field function cj> and scaled radial coordinate f. In the 
rotating case we find 

4>(f) = (j)^ 1 ((piK^ 3 <f)(r) , with f = ( (piK^ 3 r . (42) 

The ratio (w^ ax — assumes the same value ui^ for all k. For the perturbative expansion 
for the metric functions 

/ = e s ^ , m = e 5 ^ , n = e 5x ^ , u = S u w{r) , (43) 

2 

we find 6 = f^iK^J and 8 U = (piK^. Substitution in the boson field equation and the Einstein 
equations yields to lowest order in </>i 

+ »Lx) » ^ = f 3 ^™ ' M = X = , (44) 




and 



i( f3 ^)- 3 ™ =4 <-^ 2 - (45) 



Thus the limiting solution for (f> is different in the rotating and non-rotating case, as seen in 
Fig. |31 To lowest order in <jyy the mass M and charge Q reduce to the same expressions as in the 
non-rotating case, Eqs. ((36)) and (|34l) . respectively. Hence again we find 

M = uwQ = m B Q (46) 

in the limit a;. — > u m „. 



4.2.2 Stability 

Also the stability analysis for the rotating solutions is rather analogous to the case of the non- 
rotating solutions. 

When considering the mass as a function of the charge we obtain the same type of cusp 
structure, as seen in Fig.[6l A classically stable branch exists only for weak gravitational coupling. 
With increasing k, the relative extent of the classically stable branch decreases, whereas the relative 
extent of the classically unstable branch increases, until only unstable solutions remain. 

The large size of the spirals make the associated cusp structure better recognizable for the 
rotating boson stars. Still some inlets in Fig. [B] highlight this cusp structure due to the presence of 
the spirals, where at each new cusp another unstable mode arises according to catastrophe theory 

P3- 



5 Discussion 

We have addressed boson stars and their flat spacetime counterparts, Q-balls, in 5 spacetime di- 
mensions. Whereas non-rotating solutions are spherically symmetric, rotating solutions generically 
possess only axial symmetry, with two independent planes of rotation and thus two independent 
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Figure 6: Mass M versus charge Q for rotating boson stars for several values of the gravitational 
coupling constant: (a) k = 0.001, (b) « = 0.01, (c) « = 0.05, (d) k = 0.075, (e) k = 0.1; (f) 
K = 0.3; also shown is the mass of Q free bosons. 
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angular momenta. By employing a complex doublet of boson fields and choosing a special ansatz, 
where both angular momenta become equal, however, the symmetry of the solutions can be en- 
hanced. The angular dependence then becomes trivial and can be integrated, simplifying the 
problem to cohomogeneity-1. 

The rotating Q-balls and boson stars possess a quantized angular momentum. In 4 dimensions, 
for the slowest rotating states, the quantization condition reads |J| = Q. In 5 dimensions this 
condition generalizes for the equal angular momentum solutions studied to | J± \ + | J2I = 2| J\ = Q. 

As in 4 spacetime dimensions, Q-balls and boson stars exist only in a limited frequency range. 
For Q-balls we do not observe a qualitative difference in 4 and 5 dimensions. For boson stars, 
the same is true for the smaller values of the frequency, where charge and mass exhibit a spiral- 
like frequency dependence. In particular, both in 4 and 5 dimensions the effect of rotation is 
to considerably elongate the spirals for not too small gravitational coupling. However, when the 
maximal frequency aj max is approached, which corresponds to the mass of the bosons, boson stars 
show a qualitatively different limiting behaviour. 

In 4 dimensions both mass and charge tend to zero for boson stars as w s — > w max . In contrast, 
in 5 dimensions both mass and charge of boson stars tend to finite values in this limit. Employing 
scaled functions and a scaled radial coordinate we have analyzed this limit explicitly, by construct- 
ing the limiting solutions. The different behaviour can then be traced back to the different powers 
of the radial coordinate in the volume element in 4 and 5 dimensions. One may anticipate, that 
in dimensions greater than 5 the further powers of the radial coordinate in the volume element 
lead to a divergent behaviour of the mass and charge. Indeed, we have verified for non-rotating 
boson stars in 6 dimensions, that 4» and v have the same scaling behaviour and thus the mass and 
charge diverge as lu s — > w max . 

The classical stability of Q-balls and boson stars can be analyzed according to catastrophe 
theory [14] . implying a change of classical stability at each cusp encountered, when the mass 
is considered as a function of the charge. While the stability of Q-balls does not change when 
going from 4 to 5 dimensions, the stability of boson stars must be reconsidered, because the cusp 
structure changes close to w max . Thus whereas in 4 dimensions there is always a classically stable 
branch, with possible astrophysical relevance, the boson star solutions in 5 dimensions may not 
possess any classically stable branch, when the gravitational coupling becomes large. 

While we have focussed on the fundamental branches of boson stars, we have also obtained 
excited (unstable) branches, where the scalar field function has up to 6 nodes. A systematic study 
of these excited branches is still missing, however. 

Another point of interest in the 4-dimensional case is the possible occurrence of an ergoregion, 
since this would imply an instability, associated with superradiance [20] . Our previous analysis 
in 4 dimensions showed that ergoregions can be present for boson stars on the classically stable 
branch. For boson stars in 5 dimension we observe also ergoregions in some parameter ranges. 
Note that in the case of two equal angular momenta the ergoregion forms a shell. 

It is straightforward to extend the present calculations to non-rotating Q-balls and boson stars 
in more than 5 dimensions. For the rotating boson stars, however, one needs to restrict to odd 
dimensions and choose all angular momenta equal, in order to keep the problem cohomogeneity- 
1. Alternatively, allowing for general angular momenta will increase the numerical complexity 
tremendously, since with each new plane of rotation the dependence on a further angular coordinate 
arises. 
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